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A discontinuity in the low-mass IMF — the case of high 
multiplicity 
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ABSTRACT 

The empirical binary properties of brown dwarfs (BDs) differ from those of normal 
stars suggesting BDs form a separate population. Recent work by Thies and Kroupa 
revealed a discontinuity of the initial mass function (IMF) in the very-low-mass star 
regime under the assumption of a low multiplicity of BDs of about 15 per cent. How- 
ever, previous observations had suggested that the multiplicity of BDs may be sig- 
nificantly higher, up to 45 per cent. This contribution investigates the implication of 
a high BD multiplicity on the appearance of the IMF for the Orion Nebula Cluster, 
Taurus- Auriga, IC 348 and the Pleiades. We show that the discontinuity remains pro- 
nounced even if the observed MF appears to be continuous, even for a BD binary 
fraction as high as 60%. We find no evidence for a variation of the BD IMF with 
star-forming conditions. The BD IMF has a power-law index obd ~ +0.3 and about 
2 BDs form per 10 low-mass stars assuming equal-mass pairing of BDs. 

Key words: binaries: general — open clusters and associations: general — stars: 
low-mass, brown dwarfs — stars: luminosity function, mass function 



1 INTRODUCTION 

The origin of brown dwarfs (BDs) remains the subject of 
intense discussions. There are two broad ideas on their 
origi n: 1. the classical star-like formation scenario of BD s 
(e.g. |Xd ams fe Fatuzzol 1 19961 : | Padoan fc Nordlundl |2004| ). 
and 2. BDs and some very-low-mass stars (VLMSs) form 
as a separate population (hereafter named BD-like be- 
sides the classical star-like population) with a different 
format ion history than stars, e .g. as ejected stellar em - 
bryos (|Reipurth fc Clarke! I2OO1I: iKroupa fc Bouvierl l2003al) 



or as disrupted wide bi naries (jGoodwin fc Whitworthl 2007 
IStamatellos et aT]|2007l ). Additionally, the formation of BDs 
from Jeans instabilities in high-density filaments near the 
centre of a massive star-forming cloud has been recently 
suggested implying such BDs to be preferenti ally located in 
clust ers with strong gravitational potentials (jBonnell et alj 
l200Sh . 

The star-like formation scenario fails to reproduce the 
observed different binary properties of BDs and stars. Es- 
pecially the truncation of the semi-major axis distribu- 
tion between 10 and 20 AU for BDs and the different 
mass-ratio distr i bution of BDs and stars dBouv e t al. 2003; 
Burgasser et al.1 120031: iMartfn et al.1 120031 



ulation to stars. The assumption of two separate popula- 
tions would then require two separate initial mass functions 
(IMFs) of the individual bodies of a star cluster. Although 
the observed mass function may appear approximately con- 
tinuo us from the lowest mass BDs up to the highest mass 
stars (|Lodieu et al.| [2007). a discontinuity in the IMF may 
be present but be masked by 'hidden' (unresolved) bina- 
ries only emerging if the observed MF is corrected for unre- 
solved multiplicity. This issue has been discussed in greater 
detail in Thies & Kroupa (2007, hereafter TK07) for the 
case of a low multiplicity of 15 % of the BD-like population. 
However, a higher BD multiplicity betw een 20 % and 45 % 
had b e en reported by some au thors, e.g. Ijeffries fc Maxtedl 
( 2005V iBasri fc Reined | 2006). This raises the need, dealt 
with this contribution, for including higher multiplicities as 
well as for an analysis of the general effects of a high multi- 
plicity on the IMF and on the BD-to-star ratio. 



20031: 



Kroupa et al.l 
the BD desert 



.Close et al.l |2003h" as well 
I McCarthy et all2003l : lGrether fc Lineweaverl2006l ) , are dif- 
ficult to account for if BDs form indistinguishably to stars 
This implies the need of treating BDs as a separate pop- 



In Section [2] we review the evidence for a separate BD- 
like population. Section [3] briefly introduces the mathemat- 
ical method of calculating the IMF including unresolved bi- 
naries. In Section [4] the new results are presented and com- 
pared to those of TK07. The Summary follows in Section [5] 
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Figure 1. Mass ratio distribution of binaries in a Monte Carlo 
sample cluster following the canonical stellar mass function (§|3j. 
Upper panel: Binaries formed via random pairing over the com- 
plete mass range of a single population of BDs and stars (dots). 
The companions have masses down to 0.01 Mq (thin dashed line) 
for all kinds of primaries. T his is in contradiction with observa- 
tions llLafreniere ct al. 200i|, black diamonds). Lower panel: Bi- 
naries from random pairing within separate populations of BDs 
between 0.01 and 0.075 Mq, and stars above 0.075 Mq, i.e. there 
are no stellar companions below 0.075 Mq (thick dashed line). 
T he resulting mass-ratio d istribution fits well to the one observed 
by Lafrcnicrc ct al. ( 2 0081 ) in the stellar regime. The regime of the 
(approximately equal-mass companion) BD-like binaries is indi- 
cated by the shaded region. 



2 BROWN DWARFS AS A SEPARATE 
POPULATION 

2.1 Motivation 

It may be argued that by treating the BDs as a separate pop- 
ulation this forces an IMF discontinuity near the stellar/sub- 
stellar mass limit by construction. Indeed, the semi-major 
axis data and binary fraction (here used as a simplifica- 
tion for the multiplic ity, neglecting multiples of higher order, 
iGoodwin et al.ll2007l ) as a function of the primary mass can 
be interpreted to be continu ous with no evidence for BDs 
being a separate population (|Burgasser et al.ll2007l ). 

Given this argumentation, it is essential to describe the 
methodology applied in our analysis: We seek one mathe- 
matical formulation which is a unification of the binary pop- 
ulation for G-, K-, M-dwarfs and VLM-stars and BDs. This 
is found to be possible for G-, K- and M-dwarfs: thus, for ex- 



ample, G-dwarfs have mostly K- and M-dwarf companions, 
and the period-distribution functions of G-, K- and M-dwarf 
primaries are indistinguishable. Further, the mass-ratio and 
period-distribution functions for G-, K- and M-dwarf pri- 
maries derive from a single birth mass-ratio and period dis- 
tribution which does not differentiat e according to the mas s 
of the primary (|Kroupa et al.ll2003l ; [Goodwin et aj]|2007h . 
One single mathematical model can therewith be written 
down which treats G-, K- and M-dwarfs on exactly the same 
footing - one can say that G-, K- and M-dwarf stars mix ac- 
cording to one rule (random pairing from the IMF at birth). 

If BDs are to be introduced into a similar mathemat- 
ical formulation which does not differentiate between BDs 
and stars, then the model fails, because it leads to (1) a too 
wide BD period-distribution function, (2) too many BD bi- 
naries, (3) far too man y stellar-BD binari e s, and (4) far too 
few st ar-star binaries 1 Kroupa et al.ll2~003l ). lLafreniere et al.l 
(2008) show that the mass-ratio of binaries depends on the 
primary-star mass in a way that results in an almost con- 
stant lower mass limit of the companions near 0.075 Mq 
for the Chamaeleon I star-forming region (see the figures 7 
and 12 in their paper). This distribution can be reproduced 
by random pairing over the stellar mass range, while global 
random pairing over BDs and stars yields a different distri- 
bution, as shown in Fig. [1] In order to avoid these failures, 
and in particular, in order to incorporate the BD desert 
('stars and BDs don't mix, while G, K and M dwarfs do 
mix') into the mathematical formulation of the population, 
it is unavoidable to invent special mathematical rules for the 
BDs. That is, stars and BDs must be described separately. 

TK07 show that this necessarily implies a discontinuity 
in the IMF, given the observational data. We emphasise that 
the observational mass distributions lead to this conclusion, 
once the correct mathematical description is incorporated 
consistently. However, TK07 assume a rather low binary 
fraction of the BD-like population of 15 %. Some papers 
iGuenther fc Wuchterll 120031 ; iKenvon et all |2005l : Ijoergensl 
120061 ) that report the discovery of close BD binaries instead 
conclude that these may imply a si gnificantly higher binary 
fraction between 20% and 45% (jjeffries fc Maxtedl [2003 : 
iBasri fc Reiners!l2006l ), the latter being si milar to that o f 
stars in dynamically evolved environments (|Kroupalll995d ). 
It is therefore useful to re-address the problem TK07 posed 
by incorporating a larger BD binary fraction into the anal- 
ysis. 



2.2 A short review of binarity analysis 

This contribution is, like TK07, part of a series on the 
theor etical interpre t ation of obs e rvational stella r cluster 
data. iKroupa et all (|l99ll . ll99St ); iKroupal (|200lh showed 
for the first time that the true individual-body stellar 
IMF is changed significantly by correcting for the bias due 

Kroupal 



to unr esolved binary stars. A detailed study in 



1995all3lbT). in IKroupa. Petr fc McCaughreanl (|l999h . and 



Kroupa. Aarseth fc Hurley I 200ll ) of the observed binary 
properties of field stars, stars in the Orion Nebula Clus- 
ter (ONC), the Pleiades and the Taurus- Auriga association 
(TA) led to a thorough understanding of the energy distri- 
bution of binary systems. This work established that simply 
taking observed distributions can lead to wrong interpreta- 
tions, unless the counting biases and stellar-dynamical evo- 
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lution is treated systematically and consistently; the basis of 
the argument being that N ewton's laws of motion cannot be 
ignored. In a recent paper llteipurth et al.l (l2007h hav e sup- 
ported the predictions made by iKroupa et alj (|1999) con- 
cerning the binary population in the ONC, by uncovering 
a radially dependent binary fraction in nice agreement with 
the theoretically expected behaviour. The late-type stellar 
binary population is therewith quite well understood, over 
the mass range between about 0.2 Mq and 1.2 Mq. The 
above work has also established the necessity to correctly 
dynamically model observed data in order to arrive at a con- 
sistent understanding of the physically relevant distribution 
functions. 

Brown dwarfs (BDs), which extend the mass scale down 
to 0.01 M0, have been added i nto the theoretical ana l- 
ysis in IKroupa fc Bouvier! (|2003bl lah: IKroupa et all (I2003T) 
This theoretical study of observational data (IClose et al.l 
120031 ; iBouv et all l2003h showed that BDs cannot be un- 
derstood to be an extension of the stellar mass regime (as 
is often but wrongly stated). The hypothesis of doing so 
leads to incompatible statistics on the star-star, star-BD 
and BD-BD binary fractions, and on their energy distri- 
butions. This work showed that BDs must be viewed as 
a separate population, a nd the theoretical suggestion by 
iReipurth fc Clarke! l|200lf) . that BDs are ejected stellar em- 
bryos, is one likely explanation for this. In fact, their propo- 
sition logically implies different binary properties between 
stars and BDs, because ejected objects cannot have the same 
bi nding energies as not-eje c ted ob jects. Likewise, the model 
of iGoodwin fc Whitworthl (120071 ). according to which BDs 
are born in the outer regions of massive accretion disks, im- 
plies them to have different pairing rules than stars. 

TK07 and this contribution are a logical extension of 
the above findings. Here we repeat parts of the analysis of 
TK07 with assumed BD-like binary fractions up to 60 % 
as an upper limit. T he clu sters we analys e are the ONC 
ilMuench et al.l I2002T). TA jLuhmanl |2004| ; iLuhman et~aH 
I2003T I, IC 3 48 (iLuhman et alj|2003l) and the Ple i ades b ased 
on data by iDobbie et al.l (|2002T ). iMoraux et al.1 J2003I ') and 
the Prosser and Stauffer Open Cluster Databasqj. The aim 
is to check whether our previous results are still valid for a 
higher binary fraction, and how robust they are for different 
accounting of unresolved binary masses. 



3 IMF BASICS AND COMPUTATIONAL 
METHOD 

The IMFs are const ructed from pow er-law functions similar 
to that proposed by Sal peterl (|l955f ). 



£ (m) = - — — km a , 
dm 

or in bi-logarithmic form 



£ L (log 10 m) = 



dn 



dlog 10 m 



= In (10) m £(m) = k^m 



(1) 



(2) 



where fc is a normalisation constant and fcL = ln(10)fc. While 
Salpeter found a ~ 2.35, the canonical stell ar IMF, f s t ar , i s 
constructed as a two-part-power law after iKroupal (2001), 



with ai = 1.3 for a stellar mass m < 0.5 Mq and a = 2.3 
for higher masses. The substellar IMF, £bd, is taken to be 
a single power-law with cluster-dependent exponent obd- 

The basic assumption is that a large fraction of bina- 
ries remains unresolved since cluster surveys are often per- 
formed with wide-field surveys with limited resolution. One 
may be tempted to use the observed IMF (hereafter IMF b s ) 
as a direct representation of the true IMF of individual bod- 
ies (simply the IMF hereafter). However, the observed IMF 
(IMF bs) can differ largely from the IMF, especially at the 
low-mass end of the population which contains most of the 
stellar companions. If the companion has a much lower mass 
than its primary, then its light does not contribute much to 
the combined luminosity and spectral type, and thus the 
derived mass is essentially that of the primary. If, however, 
both components have near-equal masses (as expected from 
random pairing for very low primary masses), the low-mass 
region of the IMF D b s may be depressed even further, since 
the combined luminosity can be up to twice the luminosity 
of the primary alone. Therefore a fraction of unresolved low- 
mass binaries is counted as single stars, maybe even of higher 
mass, while their companions are omitted, attenuating the 
IMF b s at the lower mass end. 

Possible approximations to the IMF b s are the system 
IMF (IMFsy S ), that is the IMF as a function of system mass 
(see equations 6 to 8 in TK07), and the primary body IMF 
(IMFp r i), the IMF as a function of the primary object mass, 



£pri(f7lprim) — /tot-^bod 



£(m prim )£(m) dm, 



(3) 



where iVbod is the total number of objects, m m i n is the min- 
imum mass of an individual body in the given population, 
/tot = A r bny/(A r sng + iVbny) is the total binary fraction, and 
£(m) = £(m) / Nbod is the normalised individual-body IMF. 

In TK07 the IMF sys has been used for the fitting pro- 
cess. However, one may argue that the mass derived from the 
system luminosity is closer to the mass of the primary star 
since the luminosity is mainly given by the primary object 
and the spectral features of the companion are outshone by 
those of the primary. Therefore, the IMF pr i has been used 
as the workhorse in the current contribution. 

To obtain the true IMF from an observed mass distri- 
bution a binary correction has to be applied to each na- 
tive population (i.e. a population of objects that share the 
same formation history) the cluster consists of. This is done 
here via the semi-analytical backward-calculation method 
and x 2 minimisation against the observational data intro- 
duced in TK07. It assumes two native populations with dif- 
ferent IMFs, different overall binary fractions and different 
mass-ratio distributions (namely the two extreme cases of 
random pairing and equal-mass pairing for BDs while ran- 
dom pairing is always applied to stars). For each cluster the 
BD IMF slopes, the population ratio, 



iV s1 



(4) 



1 Available at |http://www. c fa.harvard.edu/~stauffer/opencl/ 



and the upper mass limit of the BD-like IMF, m maXj BD, arc 
to be fitted, while the lower mass limit of BDs (0.01 Mq) 
and of the star-like population (0.07 Mq) is kept constant. 
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Here, the number of BD-like and star-like objects is given 
by 



AW = / £star(m)dm, 



(5) 



respectively. The IMFs are then transformed into separate 
primary mass functions. Before being compared to the ob- 
servational MFs the fitted IMF pr i has been smoothed by a 
Gaussian convolution along the mass axis in order to sim- 
ulate the error of the mass determination (see TK07 for a 
more detailed description). This process is repeated itera- 
tively until \ 2 reaches a minimum. 

For the Pleiades the BD data do not constrain the 
power-law index, so fixed power-laws with qbd = 0.3 (the 
canonical value) and q:bd = 1 have been used here. It should 
be noted that the power-law indices are in ro ugh agreement 
with t he po wer-law index a = 0.6 deduced bv lBouvier et al.1 
(|l998h and iMoraux et alj l^OOSt ). Since they use BDs and 
low-mass stars up to 0.48 Mq while only BDs and VLMSs 
are used in our contribution these values have to be com- 
pared with caution. 

The crucial point in performing the binary correction is 
that the assumed number and mass range of a native pop- 
ulation largely affects the resulting IMF pr i. If, for example, 
only one overall population is assumed (as in the traditional 
star-like scenario for BDs and stars) but there are actually 
two separate BD-like and star-like populations with different 
mass ranges, then the binarity is corrected for wrongly at the 
lower mass end of the star-like population since a mixing of 
binary components between BDs and stars is assumed that 
does not exist in reality. Reversely, the observed (primary) 
IMF of a cluster may appear as being continuous while actu- 
ally consisting of two populations, because the discontinuity 
is masked by the interference of different probability den- 
sities of the populations in the transition or overlap region 
on the one hand and different binary fractions on the other, 
as well as being smeared out by measurement uncertainties. 
Thus, an apparently continuous IMF pr i or IMF sys may be 
related to a discontinuous IMF or, more precisely, a compos- 
ite IMF which can only be revealed by reducing the fraction 
of unresolved binaries to insignificance by high-resolution 
observations. 

The magnitude of the discontinuity, measured as the 
number ratio of BD-like to star-like objects at the hydrogen- 
burning mass limit (HBL), 7£hbl, is given by 



N BD (m « 0.075 M ) 

TZ-hbl = 7 r 

iVstar (m ~ 0.075 Mq) 



(6) 



If there is no overlap of the fitted BD-like and the star- 
like population (which is actually the case for the ONC and 
the Pleiades), TvLhbl is calculated from extrapolation of the 
BD-like IMF to the HBL. Since 7?.hbl depends on the binary 
fraction among each population, the binary fraction is varied 
from /bd = to /bd = 0.6 in order to include even the most 
extreme BD binary fraction. The unresolved stellar binary 
fraction, /star, is set to 0.4 for the ONC, IC 348 and the 
Pleiades while that of TA is assumed to be 0.8, as in TK07. 
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Figure 2. The best-fit IMFs for the ONC llMuench et alj|2002l. 

histogram) for /bd = 0.15, as displayed i n TK07, and for 
/bp = 0.45, the upper limit deduced by Ijeffries fc Maxtedl 
(2005). The primary IMFs (solid curves) are being derived from 
separated BD-like (dashed lines) and star-like (dotted lines) 
populations consis t ent w ith the ejected-embryo hypothesis of 
iReipurth fc Clarke! (l200lh or th e disk-fragmentation hypothesis 
of lGoodwin fc Whitworthj <2007ft . The assumed fraction of unre- 
solved star-like binaries is / 3 t a r = 0.4. Despite the high binary 
fraction in the lower panel the discontinuity near the stellar- 
substcllar border is only slightly reduced. The unequal binary 
fractions for different masses and populations mask the appar- 
ent discontinuities of the IMFs in the VLMS region in IMF pr i 
(solid curves). The thi n dashed histo g ram r efers to a substcllar 
peak in the data from iMuench et al.l ll2002h which may be due 
to non-physical artef acts in the mass-lum inosity relation used for 
the mass calculation llLada fc Ladall2003ft . 



4 RESULTS 

4.1 IMF fitting parameters for different BD 
binary fractions 

For illustration, Fig. f2] shows the fitted BD-like and star-like 
IMFs, (bd and £ stai and the resulting IMF pri for the ONC 
for an assumed /bd = 0.15 (upper panel) and /bd = 0.45 
(lower panel), both using equal-mass pairing for BD-like bi- 
naries and random pairing for star-like ones. Random pair- 
ing means, in this context, pairing two stars selected by 
chance from the IMF. The discontinuity (eq.[|J} between the 
BD-like and the star-like IMF becomes slightly smaller for 
higher binary fractions while the BD-like IMF slope remains 
almost constant. The discontinuity between both popula- 
tions is, however, still present. 

The top panel of Fig. [3] shows the dependency of qbd 
on /bd for the ONC, TA and IC 348 (the clusters for which 
«bd has actually been calculated from \ 2 minimisation) for 
both equal-mass pairing and random pairing. The IMFs have 
been fitted via the IMF pr i. It should be noted that obd > 1 
for /bd ^0.4 for random-pairing in IC 348, i.e. the turnover 
of the (bi- logarithmic) IMF in the substellar region vanishes, 
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Figure 3. Upper panel: The best-fit BD IMF power-law in- 
dices for the ONC (solid line), TA (dashed line), and IC 348 
(dash-dotted line) with equal-mass pairing and random pairing 
of BD binaries (dotted, narrow-dotted and double-dotted curves, 
respectively), as a function of the assumed BD binary fraction, 
/bDi fitted via the primary-object mass function for both BDs and 
stars (see text). The upper/lower uncertainty limits of «bd are 
about +0.3/ -0.3 for the ONC, +0.5/ -2.8 for TA and +2/ -0.6 
for IC 348. While «bd remains approximately constant for equal- 
mass pairing, there is a strong increase with increasing /bd for the 
random pairing case. Lower panel: The BD IMF power-law in- 
dices for the case of equal-mass pairing, but this time fitted via the 
BD system MF (and the stellar primary MF). The upper/lower 
limits of a B D arc about +0.3/ - 0.3 for the ONC, +0.7/ - 3.5 
for TA and +2.5/ — 0.6 for IC 348. In contrast to the equal-mass 
case in the upper panel, ojbd is now increasing with /bd- For 
comparison, the fixed «bd = 0.3 (short-dashed horizontal line) 
and ctBD = 1-0 (long-dashed horizontal line) for the Pleiades have 
been added. 



found to be around 0.1 Mq and between 0.15 and 0.23 Mq, 
respectively. Furthermore, our results for the best-fit 1Z pop , 
the population ratio, and the magnitude of the discontinu- 
ity can be summarised as follows: For /bd = the best-fit 
lZ pop is about 0.07 for the ONC and the Pleiades while it is 
about 0.15 for TA and IC 348. It increases for larger /bd, 
reaching about twice these values for the extreme binarity 
of /bd = 0.6. That is, if a realistic value of /bd ~ 0.2 is 
assumed, we expect about 1 BD-like body per 10 star-like 
ones for the ONC and the Pleiades, and about 1 BD-like 
body per 5 star-like bo dies for the others. T his result is re- 
markable given that e.g. ISlesnick et al.l l|2004l ) state a higher 
BD-to-star ratio for the ONC than for TA and IC 348. The 
result can be interpreted as a consequence of the large mass 
overlap of the BD-like and the star-like regime in IC 348 (and 
a moderate overlap in TA), i.e. that many BD-like bodies 
are in actually very-low mass stars and thus are counted as 
normal stars. Another issue is whether the substellar peak 
in the ONC MF (see the thin das hed histogram in Fig. [2]| 
is an artefact (|Lada fc La da 2003) or a real feature. In the 
latter case, !Z pop would be significantly higher (about 75 %, 
given the histogram data) for the ONC than suggested by 
our results. 

One may criticise the way of assigning a mass to an 
observed system. In TK07 the model-observed IMF has 
been created by simply adding the masses of all compo- 
nents, i.e. IMF b s = IMF sys . Because the observed data 
are being derived from luminosity functions rather than 
from mass functions, the correct way would be to con- 
vert luminosities into masses via the mass-luminosity rela- 
tion (MLR). This would require rather complicated calcula- 
tions because full-scale modelling would involve age-spreads 
and age-dependent MLRs with very significant uncertainties 
jWuchterl fc Tscharnuterll2003r i. 

Instead, a simpler way to at least embrace the real re- 
lation is to repeat the analysis (or parts of it) by using the 
primary mass instead of the system mass. This corresponds 
to the extreme case that the contribution of less-massive 
companions is negligible. This method has been used in the 
present contribution with similar results as in TK07. In ad- 
dition, similar calculations have been made for a substellar 
system IMF and for a stellar primary IMF, for equal-mass 
pairing of BDs, and for random pairing of stars. 



although the discontinuity remains (see Section [4.2[) . Simi- 
larly, the lower panel shows the trends with /bd in the case 
of equal-mass pairing if the BD-like IMF sys is used for fit- 
ting. 

The most remarkable feature is that ckbd remains al- 
most constant for equal-mass pairing in BD-like binaries. For 
random-pairing obd increases with /bd in a similar way for 
all three clusters. A similar growth is found even for equal- 
mass pairing if IMF sys is used for fitting. For comparison, 
the constant values assumed for the Pleiades are shown in 
the lower panel of Fig. [3] (straight dashed lines at qbd = 0.3 
and «bd = TO). 

The fitting of m ma x,BD yields values slightly below 
0.07 M Q for the ONC and the Pleiades. This is probably 
due to the Gaussian smearing of log m that has been used for 
smoothing the fit. For TA and IC 348, however, m maXi BD is 



4.2 The discontinuity in the low-mass IMF 

A measure for the discontinuity at the HBL, 7?,hbl, is given 
by eq. [6] For a continuous IMF 7£hbl = 1, while values sig- 
nificantly different from 1 indicate a discontinuity. Fig.|4]dis- 
plays 7?.hbl as a function of /bd ■ For all clusters 7?.hbl shows 
a similar steady increase. The uncertainties (not shown in 
the graph) can be estimated from those of lvbd and are about 
±40 % for each value. Thus the discontinuity between the 
BD-like and the star-like IMF persists for both the system 
and the primary-body IMF for the BD-like population. It is 
largest (i.e. TZhbl is smallest) for /bd = 0. The results from 
TK07 show that this even holds true if the system IMF fit 
is applied to the stellar population. 
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Figure 4. The ratio of BD-like to star-like bodies at the 
hydrogen-burning mass limit (HBL) as a function of the BD- 
like binarity, /bd, for the ONC (solid curve), TA (dashed curve), 
IC 348 (dash-dotted curve), as well as for the Pleiades for obd = 
0.3 (thin dashed curve) and ctbd = 10 (thin dotted curve). For all 
clusters there is a similar trend towards a higher 7?-hbl with in- 
creasing binary fraction. But even for the highest plausible binary 
fraction 7?-hbl < 0.5. A continuous IMF would require 7^hbl = 1- 




log 10 n/pc" 

Figure 5. The power-law index obd and the BD-to-star ratio ft 
from our model (assuming equal-mass pairing for BDs) and from 
the literature are plotted in dependence of the logarithmic central 
ste llar density of T A, IC 3 48 and the ONC . The references for ft 
arelLuhmanl 120061) fo r TA. lPreibisch et alj ll2003h for IC 348, and 
ISlesnick et aiT l2004) for the ONC. 



4.3 IMF slope and BD-to-star ratio in relation to 
the stellar density 

Certain theories of BD formation (e.g. iBonnell et alj [20081 ) 
suggest a dependency of the rate of BD formation on the 
star-cluster density. Correlating the BD IMF index, ojbd, 
and the BD-to-star ratio, TZ, against the stellar density, n, 



may uncover s uch expected depen dencies. For TA, n = 1 — 10 
stars per pc 3 dMarti'n et al.ll200ll ). n « 500 pc" 3 for IC 348 
dDucheneet al.lll999l). and n s s 20000 pc~ 3 for the ONC 
l|Hillenbrand fc Hartmannlll998l ). The values of «bd found 
in this study are plotted against log 10 n in Fig. [5] (upper 
panel). In addition, the BD-to-star ratio, TZ, is shown in the 
lower panel, where 



n 



N{0.02M Q <m< 0.075 Mq) 



(7) 



iV(0.15M Q < m < 1M Q ) ' 

The m ass li mits are chosen in accor dance with lKroupa et all 
(2003) and iThies fc Kroupal l|2007l ). The crosses connected 
with solid lines show the results of our modelling while 
the open circle s with dashed lines a r e valu es taken from 
Luhmanl d2006f) (TA ). [preibisch et ail (|2003l ) (IC 348), and 
Slesnick et aiT(|2004D (ONC). 

For obd a regression line has been calculated. However, 
only three clusters have been analysed in this study, and 
there are large uncertainties. Especially for TA and IC 348 
the confidence range is rather large here. Thus, the linear 
fit is only poorly constrained and well in agreement with a 
constant qbd- Furthermore, TZ also does not show a signifi- 
cant trend with increasing stellar density. From our analysis 
(Fig. [5]l it follows that obd ~ 0.3 and TZ m 0.2 for equal- 
mass pairing of BDs. 



5 SUMMARY 

A discontinuity in the IMF near the hydrogen burning mass 
limit appears if the binary properties of BDs and VLMSs 
on the one hand, and of stars on the other, are taken into 
account carefully when inferring the true underlying single- 
object IMF. This implies that BDs and some VLMSs need 
to be viewed as arising from a somewhat different formation 
channel than the stellar formation channel, but this result 
has been obtained by TK07 under the assumption that BDs 
have a binary fraction of only 15 %. A higher binary fraction 
may close the gap between the stellar and the BD IMF. 
We refer to BDs and those VLMSs formed according to the 
putative BD channel as "BD-like" bodies, whereas stars and 
those BDs formed according to the stellar channel as star- 
like. The BD-like channel remains unknown in detail, but 
theoretical ideas have emerged (Sections [T1 12.21 and !4.3p . 

Here we have extended the analysis of TK07 for BD-like 
binary fractions up to 60 % for the Orion Nebula Cluster, the 
Taurus- Auriga association, IC 348 and the Pleiades by using 
slight modifications of the techniques introduced in TK07. 

As a main result, we found that the discontinuity that 
comes about by treating BDs/ VLMSs and stars consistently 
in terms of their observed multiplicity properties remains 
even for the highest BD binary fraction. These results sug- 
gest that the BD binary fraction, /bd, is not the dominant 
origin of the discontinuity in the IMF, and that, conse- 
quently, two separate IMFs need to be introduced. 

It is re-emphasised that by seeking to mathematically 
describe the BD and stellar population in terms of the rele- 
vant mass- and binary distribution functions, it is unavoid- 
able to mathematically separate BDs and VLMSs from stars. 
The two resulting mass distributions do not join at the 
transition mass near 0.08 Mq. The physical interpretation 
of this logically stringent result is that BDs and VLMSs 
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follow a different formation history or channel than stars. 
This result is obtained independently by theoretical con- 
sideration of star- formation processes (Reipurth & Clarke 
200ll; iGoodwin fe WhitworthlboO^ ; IStamatellos et al.ll2007t 
Bonnell et al.ll2008t ). 



With this contribution we have quantified how the 
power-law index of the BD-like IMF and the BD-to-star ratio 
changes with varying binary fraction of BD-like bodies. The 
BD-like power-law index, qbd ~ 0.3, remains almost con- 
stant if equal-mass pairing of BD-like binaries is assumed, 
while obd increases somewhat with increasing /bd in the 
case of random pairing over the BD-like mass range. All val- 
ues of QfBD are between —0.1 and +1.3. We also find that 
although the stellar density differs from a few stars per pc 3 
(TA) to about 20000 stars per pc 3 , the resulting qbd is 
constant within the uncertainties. Similarly, the BD-to-star 
ratio does not show a trend with increasing stellar density. 
This suggests the star-formation and BD-formation outcome 
to be rather universal at least within the range of densities 
probed here. 
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